Bound states in the continuum in graphene quantum dot structures 
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The existence of bound states in the continuum was predicted at the dawn of quantum mechanics 
by von Neumann and Wigner. In this work we discuss the mechanism of formation of these exotic 
states and the feasibility to observe them experimentally in symmetrical heterostructures composed 
by segments of graphene ribbons with different widths forming a graphene quantum dot. We identify 
the existence of bound states in the continuum in these graphene quantum dot systems by means 
of local density of states and electronic conductance calculations. 



The new material denominated graphene is a single 
layer of carbon atoms which can be fabricated by differ- 
ent methods like mechanical peeling or epitaxial growth 
Nanoribbons are stripes of graphene which can 
be obtained through high-resolution lithography by 
controlled cutting processes @ or by unzipping multi- 
walled carbon nanotubes 0- The electronic behavior of 
all these nanostructures is mainly determined by their 
geometric confinement which allows the observation of 
quantum effects such as quantum interference effects, res- 
onant tunneling and localization effects. The possibility 
to control these quantum effects, by applying external 
perturbations to the nanostructures or by modifying the 
geometrical confinement, could be used to develop new 
technological applications, such as graphene-based com- 
posite materials j§] , molecular sensor devices 0, [l(| and 
nanotransistors 

An interesting feature exhibited by certain confined 
nanostructures, such as quantum dots systems, is the 
presence of bound states in the continuum (BICs) . Their 
existence was predicted at the dawn of quantum mechan- 
ics by von Neumann and Wigner |12j | for certain spa- 
tially oscillating attractive potentials, for a one-particle 
Schrodinger equation. Much later, Stillinger and Her- 
rick [131 ] generalized von Neumann's work by analyzing a 
two-electron problem, they found BICs were formed de- 
spite the interaction between electrons. The occurrence 
of BICs was discussed in a system of coupled Coulom- 
bic channels and, in particular, in a Hydrogen atom in a 
uniform magnetic field [IH • BICs have also shown to be 
present in the electronic transport in mesoscopic struc- 
tures fl5rl2f]j ]. More recently, exploiting the analogy be- 
tween electronics and photonics, Marinica et al. 1211 . 
Bulgakov and Sadreev [22| and Prodanovic et al. [23[ 
reported the presence of BICs in photonic systems. Sev- 
eral mechanisms of formation of BICs in open quantum 
dots (QDs) have been reported in the literature. The 
simplest one is based on the symmetry of the systems 
and, as a consequence, in the difference of parity be- 
tween the QD eingenstates and the continuum spectrum 
[j^ . Another mechanism takes into account a nonzero 
coupling between bound states in the QD and the con- 



tinuum spectra. The formation of BICs would be the 
result of a destructive interference process of these res- 
onances, for certain variations of the physical parame- 
ters of the QD 25l-27|. A third mechanism for the BICs 
formation in optics, is associated with the Fabry-Perot 
interferometer [30j • 

Until nowadays, there is only one experimental work, 
reported by Capasso and co-workers [3l|, in which BICs 
were measured in semiconductor heterostructures grown 
by molecular beam epitaxy. Thereby, the search of new 
systems which could be able to reveal the existence of 
BICs, with the possibility to be measured, is a very in- 
teresting and relevant field of research. The experimen- 
tal feasibility exhibits by graphene based systems and 
the great advances in the controlled manipulation and 
measurements reported in graphene, together with the 
possibility of modifying their electronic properties by ap- 
plying external potentials, suggests that BICs could be 
observable in graphene quantum dots heterostructures. 




FIG. 1. Schematic view of a GQD structure with leads of 
width Nl = 9, a conductor region composed by two symmet- 
rical junctions of width Nb = 21 and length Lb = 3 separated 
by a central structure of length Lc = 4 and width Nc ~ 9 

In this work we study the formation of BICs in 
quantum-dot like structures, formed by segments of 
graphene ribbons with different widths connected be- 
tween each other [33} . We identify the presence of BICs in 
these symmetrical graphene quantum dots (GQDs) and 
we discuss the mechanism for their formation. We found 
that the GQD local density of states as a function of the 
energy shows the presence of a variety of sharp peaks that 
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we demonstrate are BICs. The linear conductance also 
shows the presence of resonant states which contribute to 
the electronic transmission. By changing the geometrical 
parameters of the structure, it is possible to control the 
number and position of these resonances as a function of 
the Fermi energy. 

A schematic view of the considered systems is pre- 
sented in figure [TJ The conductor is formed by two 
symmetric crossbar junctions of widths Nb and length 
Lb, and a central region that separates the junctions, 
of width N c and length Lc- Two semi-infinite leads of 
width Nl = N c are connected to the ends of the cen- 
tral conductor. We studied the different electronic states 
manifested in the system as a function of the geometrical 
parameters of the GQD structure. 

Systems are described by using a single 7r-band tight 
binding Hamiltonian, taking into account first nearest 
neighbor interactions with a hopping parameter 70. We 
consider hydrogen passivation by setting a different hop- 
ping par ameter for the carbon dimmers at the ribbons 
edges [34], j e d ge = 1-12jq. To calculate electronic proper- 
ties of the system we adopt the surface Green's functions 
matching formalism In this scheme, we divide 

the heterostructure in three parts, two leads composed 
by semi-infinite pristine graphene nanoribbons, and the 
conductor region composed by two nanoribbon crossbar 
junctions, as it is shown in figure [TJ 

In the linear response approach, the electronic con- 
ductance is calculated by the Landauer formula. In 
terms of the conductor Green's functions, it can be writ- 
ten as 37] G = ^f(E) = ^Tr [T L G%T R G£] , where 
T (E) , is the transmission function of an electron cross- 
ing the conductor region, T L / R = i[T, L / R — ^ L / R ] is 
the coupling between the conductor and the respective 
leads, given in terms of the self-energy of each lead: 
^l/rVc,l/r9l/rV l /r,c- Here, V c , l /r are the coupling 
matrix elements and Ql/r, is the surface Green's func- 
tion of the corresponding lead [l(J. The retarded (ad- 
vanced) conductor Green's functions are determined by 
37 :(;{': 1 = [E-H C - - £f A ] -1 > where H c is the 
Hamiltonian of the conductor. 

The Fig [2] displays results of the local density os states 
(LDOS) (upper panel) and the linear conductance (lower 
panel) for a GQD structure formed by two armchair rib- 
bons leads of width Nl — 5 and a conductor region 
composed by two symmetric crossbar junctions of width 
Nb = 17, length Lb — 3 and relative distances between 
the junctions Lc — 5. The conductance of a pristine 
Nl = 5 armchair nanoribbon has been included for com- 
parison (light green dotted line). 

It can be observed in the LDOS and in the conductance 
curves a series of peaks at determined energies. This res- 
onant behavior of the electronic conductance arises from 
the interference of the electronic wave functions inside 
the structure, which travel forth and back forming sta- 
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FIG. 2. LDOS (upper panel) and conductance (lower panel) 
as a function of the Fermi energy for a GQD structure based 
on leads of width Nl = 5, two symmetric crossbar junctions 
of width Nb = 17 and Lb = 3. The central region has a 
width Nc = 5 and length Lc = 5. Marks (a), (b) and (c) 
denote position of peaks in the LDOS which are absent in the 
conductance. These states are identified as BICs. 



tionary states in the conductor region (well-like states). 




FIG. 3. Corresponding contour plots of some sharp LDOS 
resonances marked in Fig. [2] LDOS (upper panel). 

We start our analysis focusing on the sharp states 
present in the LDOS displayed in this figure. We have 
marked the first three sharp states with the letters (a), 
(b) and (c). We note that the corresponding peaks in 
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the conductance are absent and we identify these states 
as BICs. We calculated the spatial distribution of these 
states, representing by the corresponding contour plots 
exhibited in the figure[31 We observe that BICs are com- 
pletely localized at the region defined by the crossbar 
junctions. Any electron from the leads, injected at these 
energies, will be spatially confined in the junctions due 
to the practically infinite (zero resonance width) lifetime 
of these states. Therefore BICs do not play any role in 
the transport properties of these GQD structures. This 
can be seen in Fig. [5J where the corresponding resonant 
peaks are absent in the conductance curves. 

We note that the bound states in the continuum ex- 
hibit a spatial symmetry in the transverse direction of 
the conductor, with the presence of nodes and maxi- 
mum in the amplitude of the LDOS along that direc- 
tion. This behavior is a consequence of the crossbar 
junction symmetry with respect to the longitudinal axis 
of the GQD structure. The localized states belonging 
to the crossbar region could interact between each other 
through the continuum of states of the pristine ribbon 
leading to the formation of bound-states in the contin- 
uum. In this sense, one of the mechanisms of forma- 
tion of BICS in our systems correspond to the first one 
described in the introduction of this letter. Follows the 



analysis of Moiseyev 32] , the number of BICs can be con- 
trolled by varying the gate potential applying in the up 
and down barriers of the GQDs. In order to get a better 
understanding of this phenomenon, we introduce a sim- 
ple model that captures the essence of the formation of 
BICs in our GQD structures. The model consists in two 
side-coupled impurities attached to a perfect quantum 
wire [1^] as shown in Fig. [4] 



density of states of the entire system is obtained summing 
over a and can be written approximately as, 
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7T (w - £ ) 2 + 4 7 2 7T (W - £ ) 2 + (<5 2 /2 7 ) 2 



(1) 



The density of states is then the sum of two Lorentzian 
shapes lines with widths T+ = 2 7 and T_ = S 2 /2^, cor- 
responding to those states strongly and weakly coupled 
to the continuum, respectively. In the limit of S — > 0, 
r_ vanishes and the state weakly coupled to the con- 
tinuum becomes a bound-state in the continuum. This 
state arises from the interference of the localized states 
in the impurities, through the continuum states of the 
wire. In this sense, this interference phenomenon is simi- 
lar to phenomena like Fano and Aharonov-Bohm effects. 
If we denote as tp u and tpd the wave functions of the up 
and down impurity respectively, it is straightforward to 
show that the antisymmetric state, (ip u — ipd) is an eigen- 
state of the complete system and therefore it is a bound- 
state in the continuum. In analogy, in the GQD structure 
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FIG. 4. Scheme of side-coupled impurities attached to a per- 
fect wire. 

By using the Dyson equation G = g+gVG we calculate 
the Green's function (G) in terms of the corresponding 
Green's function of the isolate sub-systems (g), here V is 
the matrix coupling between the impurities and the wire 
(Vqu = Vod = Vo). To obtain the LDOS of each impurity, 
Pa , (a — u,d), we calculate the imaginary part of the 
diagonal elements of the Green's functions, G a - Setting 
the site energies as, e u = £0 + $ and £d = £0 — S, 7 = 
TrV 2 p(0), where p(Q) corresponds to the LDOS in the site 
of the wire without impurities, and taking S <C 7, the 



FIG. 5. LDOS (upper panel) and conductance (lower panel) 
as a function of the Fermi energy for a GQD structure com- 
posed by the same parameters of Fig. [2l with an up-down 
asymmetric gate potential applied to the junction regions. 

the formation of the BICs follows the same mechanism. 
According to it, if any infinitesimal small perturbation 
breaks the transversal symmetry, the BICs become res- 
onant states with infinitesimal widths. The widths of 
these "quasi-BICs" can be controlled, for example, by 
tuning the asymmetry of the system through gate po- 
tentials. For instance, figure [5] displays the LDOS and 
conductance as a function of the Fermi energy, for the 
same systems considered in figure [5J but now taking into 
account a small up-down asymmetric gate potential ap- 
plied to the edges of the GQD. Due to this external per- 
turbation, it is possible to observe the apparition of new 
peaks of conductance at the BICs energies levels in com- 
parison with the unperturbed case. Therefore, for the 
feasibility to observe this phenomenon, it is necessary to 
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build a GQD as symmetrical as possible and to control 
the asymmetry via gate potentials. By measuring the 
changes in the widths of the quasi-BICs in the conduc- 
tance, it would be possible to obtain indirect evidence of 
the BICs. 

In summary, in this work we have studied the forma- 
tion of bound-sates in the continuum in quantum-dot like 
structures. We identify the presence of these states in the 
LDOS in symmetrical graphene quantum dots structures 
and we discuss the mechanism of formation of these kind 
of exotic states. Our results suggests that BICs could be 
observable in GQDs. The BICs or quasi BICs can have 
applications in the fields of the spintronics and the quan- 
tum computing. For instance, by applying a magnetic 
field to the GQD structure, the spin degeneracy is lifted 
and by controlling the resonances width with a gate po- 
tential, an efficient spin filter can be designed. On the 
other hand, the presence of two simultaneous BICs could 
be used as a qubits for quantum information. Also, BICs 
could have important applications in photonic crystals, 
so in this sense, we expect that the study of this kind of 
states in graphene nanoribbons quantum dots could open 
a new line of applied research. 
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